Multiplicity Fluctuations and Correlations in 
Limited Momentum Space Bins in Relativistic Gases 

Michael Hauer/'^ Giorgio Torrieri,^ and Spencer Wheaton^ 

^ Helmholtz Research School, University of Frankfurt, Frankfurt, Germany 
^ UCT-CERN Research Centre and Department of Physics, 
University of Cape Town, Rondebosch 7701, South Africa 
^ . ^FIAS, J. W. Coethe Universitdt, Frankfurt am Main, Germany 

^ ■ Abstract 

■ Multiplicity fluctuations and correlations are calculated within thermalized relativistic ideal quantum 
^ ', gases. These are shown to be sensitive to the choice of statistical ensemble as well as to the choice of 
^SJ ' acceptance window in momentum space. It is furthermore shown that global conservation laws introduce 
non-trivial correlations between disconnected regions in momentum space, even in the absence of any dy- 
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I. INTRODUCTION 



Fluctuations of, and correlations between, various experimental observables are believed to have 
the potential to reveal new physics. The growing interest in event-by-event fluctuations in strong 
interactions is motivated by expected anomalies in the vicinity of the onset of deconfinement 
a„d in the case when the^expanding system goes through the transition hne between quaric-gluon 
plasma and hadron gas [5]. In particular, a critical point of strongly interacting matter may be 
accompanied by a characteristic power-law pattern in fluctuations GJ. Recently, it has been suggested 
that correlations across a large interval of rapidity arise from a color glass condensate M, Nl . In recent 

n n 

years a wide range of experimental measurements of fluctuations of particle multiplicities [9|, |10| . 




have been reported. 



transverse momenta and multiplicity correlations in rapidity [12 
leading to a lively discussion regarding their physical interpretation [8 

To get a reliable indication of new physics, it is important to note that most fluctuation and 
correlation observables are also sensitive to some "baseline" contributions that, nevertheless, can 
have non-trivial behaviour. For instance, most fluctuation and correlation observables are sensitive 
to the global characteristics (e.g. the distribution of the number of colliding nucleons) of a sample of 
events, which can in turn be non-trivially constrained by centrality bin construction [15]. Similarly, 
conservation laws can provide a "trivial" correlation between observables. The effects of such corre- 
lations depend on the scale at which these conservation laws become important. This scale could be 
an5d;hing, from microscopic (mean free path, diffusion scale) to the macroscopic size of the system. 

The purpose of this paper is to study these baseline correlations in a limiting case: that of 
a thermalized relativistic ideal (no inter-particle interactions) quantum gas for which we want to 
assess the importance of globally applied conservation laws for particle multiplicity fluctuations and 
correlations. In this case, all observables are calculable simply using statistical mechanics techniques. 
Such an approach has a long and distinguished history of calculating particle multiplicities in hadronic 



collisions 
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Conventionally in statistical mechanics three standard ensembles are discussed; the micro canon- 
ical ensemble (MCE), the canonical ensemble (CE), and the grand canonical ensemble (GCE). In 
the MCE^ one considers an ensemble of micro states with exactly fixed values of extensive conserved 
quantities (ener gy, momentum, electric charge, etc.), with 'a priori equal probabilities' of all micro 



states (see e.g. 



301]). The CE introduces the concept of temperature by introduction of an infinite 



thermal bath, which can exchange energy (and momentum) with the system. The GCE introduces 
further chemical potentials by attaching the system under consideration to an infinite charge bath^. 



^ The term MCE is also often applied to ensembles with energy but not momentum conservation. 

^ Note that a system with many charges can have some charges described via the CE and others via the GCE. 



Only if the experimentally accessible system is only a small fraction of the total, and all parts have 
had the opportunity to mutually equilibrate, can the appropriate ensemble be the Grand Canonical 
one. 

In the limit of very large volume and constant density (the thermodynamic limit), average values 
of intensive quantities are the same for all ensembles. However, even in this limit, these ensembles 

n 

have different properties with respect to fluctuations and correlations [31|- In the MCE, energy and 
charge are exactly fixed. In the CE, charge remains fixed, while energy is allowed to fluctuate about 
some average value. Finally in the GCE the requirement of exact charge conservation is dropped, too. 



One may also consider isobaric ensembles 
33I . l34| . In previous articles 
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or even more_general 'extended Gaussian ensembles' 
3^ 



35 



36 



33, 



39 



40 



411 ] it was shown that these 



differences mean that multiplicity fluctuations are ultimately ensemble specific. 

In this article we extend these results to fluctuations and correlations between particle multiplic- 
ities in limited bins of momentum space (rapidity transverse momentum pt and azimuthal angle 
(f)). In section [IT] we present details of the calculation of correlations within statistical mechanics. The 
following two sections present calculated fluctuations and correlations within the same momentum 
space bin for a stationary (section [III Al) and boosted (section IIIIBI) system. In section HVl we discuss 
long range correlations between momentum space bins. A discussion section |Vl summarising our 
results and discussing their phenomenological implications within the context of heavy ion collisions, 
closes the paper. 



II. CORRELATIONS AND FLUCTUATIONS WITHIN DIFFERENT ENSEMBLES 



In a recent paper |4l| we have shown that GCE joint distributions of extensive quantities con- 
verge to Multivariate Normal Distributions (MND) in the thermodynamic limit (TL). MCE or CE 
multiplicity distributions could then be defined through conditional GCE distributions. In general 
one may write for the multiplicity distribution Pce{N) of a CE with conserved electric charge Q: 

p number of all states with Q and 

'^'^ number of all states with Q 

Likewise one can write for the CE joint multiplicity distribution of particle species A and B: 

P (N N ) iiumber of all states with Q, Na and Nb 
' number of all states with Q 

The number of all micro states with electric charge Q, and multiplicities A^^ and Nb of a system 
with temperature T and volume V is given by the CE partition function Zce{V,T,Q, Na, Nb)- 
Similarly, Zce{V, T, Q) denotes the number of micro states with fixed electric charge Q, but arbitrary 
multiphcities A^i and Nb, for the same physical system. 



The strategy to calculate joint multiplicity distributions could thus be the following (in principle 
also valid at finite volume): 

p ^ . _ ZMV,T,Q,N^,Ns) 

P-^^^^^'^ - Zce{V,T^Q) ' 

^ e^T Zce{V,T,Q,Na,Nb) ZGCEiV,T,fx) 

ZGCE{V,T,fi) e^^ Zce{V,T,Q) ' 

= PgceiQ,NA,NB) P-i(Q) = PgceiNA,NB\Q) . (5) 

In order to get from Eq.([n]) to Eq.([5]) both canonical partition functions ZcE(y,T,Q, Na, Nb) and 
Zce{V,T,Q) are divided by their GCE counterpart ZccEiY^T, ^) and multiplied by e'^^ . The 
first term on the right hand side of Eq.(jl]) then equals the GCE joint distribution PgceiQ, Na, Nb), 
while the second term is just the inverse of the GCE charge distribution PgceiQ)- Their ratio is 
the (normalised) GCE conditional distribution of particle multiplicities Na and Nb at fixed electric 
charge Q, Pgce{NA, Nb\Q), and equals the CE distribution Pce{NA, Nb) at the same value of Q. This 
result is independent of the choice of chemical potential fi. 

The problem of finding a solution, or a (large volume) approximation, to the CE distribution 
Pcei^A, Nb) is now turned into the problem of finding a solution or approximation to the GCE 
distribution of multiplicities Na and Nb, and charge Q. The role of chemical potential (or Lagrange 
multiplier) /i will be discussed in Section IIIII 

From the assumption that the GCE distribution PgcdQ, Na, Nb) converges to a Tri-variate Nor- 
mal Distribution, it also follows that the marginal distribution PgceiQ), well as the conditional 
distribution Pgce{NA, Nb\Q), are Normal Distributions. Hence, Pce^Na, Nb) should have a good ap- 
proximation in a Bivariate Normal Distribution (BND) in the large volume limit (where the large 
particle numbers can be appropriately treated as continuous): 

Pbnd{Na,Nb) = , \ (6) 

2nV^/aiaj^ (1 - p2) 



X exp 



1 

2V 



(ANa)^ ^ (ANa) (ANb) ^ (ANb)' 



crlil-p') "^aAaBil-p") ^1(1 - p') _ 
where ANx = Nx - (Nx), with X = A,B and: 

Val ^ (iVl) - {Na)' , (7) 
Val ^ (iV|) - {Nb)' , (8) 
VaAB = {NaNb) - {Na){Nb) . (9) 

Here VaA and VaB are the variances of the marginal multiplicity distributions of particles Na and 
Nb- The term Vctab is called the co- variance. Additionally we define the scaled variance ujx'- 

{Nl) - {Nx)' , , 



which measures the width of the marginal distribution P{Nx). Lastly, 

P = , (11) 

O'aO'B 

is the correlation coefficient between particle numbers Na and Nb- 

The distribution Eq.Q hence has 5 parameters: Mean values (Na) and (Nb), variances of 
marginal distributions Va\ and Vcr^, and the correlation coefficient p. Loosely speaking, the corre- 
lation coefficient p defines how the BND, Eq.([6]), is tilted. In the case where p > 0, the distribution is 
elongated along the main diagonal, and measuring a larger (smaller) number of particles Na implies 
that it is also more likely to measure a larger (smaller) number of particles Nb- The distribution 
is tilted the other way, if p < 0. In this case, multiplicities Na and Nb are anti-correlated, and 
measuring A^^ > (Na) implies that it is now more likely to measure Nb < (Nb). Particle numbers 
Na and Nb are uncorrelated, if p = 0. 

Similarly, we define MCE multiplicity distributions in terms of conditional GCE distributions of 
extensive quantities. For this we will ffist find a suitable approximation to the GCE joint distribution 
of extensive quantities (electric charge, energy, momentum, particle number(s), etc.) by Fourier 
analysis of the GCE partition function. The MCE multiplicity distribution is then given by a slice 
along a surface of constant values of extensive quantities. 



A. GCE Partition Function 



The GCE partition function of a relativistic gas with volume V, local temperature T = chem- 
ical potentials p and collective four velocity u reads (the system four-temperature [42] is P = [3u): 



'GCE 



(F, (3, p, u) 



exp 



(12) 



where \& (/3, p, u) is a sum over the single particle partition functions ipi iP, P, u) of all particle species 
considered in the model: 

^{P,p,u) = V'^(/3,P,w) . (13) 



The single particle partition function ipi (/3, p, u) of particle species 'Z' is given by a Jiittner distribu- 
tion: 



V'z {P, P, u) 



91 



j £p ln(^l ± e-f^P" + 



±1 



(14) 



where pf are the components of the four momentum, qj are the components of the charge vector and 
gi is the degeneracy factor. The upper sign refers to Fermi-Dirac statistics (FD), while the lower sign 
refers to Bose-Einstein statistics (BE). The case of Maxwell-Boltzmann (MB) statistics is analogous. 



In the following we restrict ourselves to systems moving along the z-axis and use variables pt 
and (f). For a boost in rapidity of one finds for the four-velocity u, the four-momentum pi and the 
integral measure d^p, respectively: 



u = (cosh?/o, 0, 0, sinhyo) , (15) 
Pi = ^y^m^~T~p|^cosh?/, Pt cos (p, pTsincj), y^?T7|~+~p|^ sinh , (16) 



d^p = Pt +Pt cosh(?/ — y^) dy dpx d(f) , (17) 

where mi is the mass of a particle of species /. The single particle partition function Eq. flT^ now 
reads: 

7r OO OO 

ilJi{f3,fl,u) = jdcpj dpT j dy pTsJm1 + plcosh{y-y,) In (l ± g"^ «m + m.' 



-TT — OO 

(1^ 



For the examples in the following sections we chose a simple gas with only one conserved charge, 
denoted as a 'pion gas'. The presented formulae are, however, also readily applicable to a hadron 
resonance gas (HRG). Depending on what system one may want to study, one introduces chemical 
potentials /I and the 'charge' vector of particle species /: 

jl = {fiB, fJ'S, fJ'Q, fJ'NA, f^Ng) qi = (bi, si, qi, nA{^), riBi^l)) for a HRG , (19) 
fl = (/iQ, fiNA, f^Ng) qi = {qi, nA{fl), ^^(fi)) for a pion gas , (20) 

where fiB, fJ'S, and /ig are the baryon, strangeness, and electric charge chemical potentials, respec- 
tively. fiN^ and /iiVfl are particle-specific chemical potentials, and could denote out of chemical 
equilibrium multiplicities of species 'A' and 'B', similar to phase space occupancy factors 75 4J] and 
7g |45|. Throughout this paper we neglect finite density effects, so = = 0. 

In addition, bi, si, and qi are the baryonic charge, the strangeness, and the electric charge of 
particle species '/'. Q is the momentum space bin in which we wish to measure particle multiplicity. 
nyi(i7) = 1 if the momentum vector of the particle is within the acceptance, n^(fi) = if not. The 
charge vector g] also contains, to maintain a common notation for all particle species considered in 
Eq.(IT3!l. the 'quantum' number Ub^^)- 

One may also be interested in correlations of, for instance, baryon number B and strangeness S, 



as e.g. in Refs.[46|, |47|. In this case, the A particle, with = (1, —1,0, 1, 1), would be counted in 



groups A and B, provided the momentum vector is within the acceptance Q. 



B. Generating Function 



To introduce the generating function of the charge distribution ZccEiV, P, u; 0, a) in the GCE 
we substitute in Eq.([Ti 



3 ' 



13 ^ j3 



la 



(21) 
(22) 



The yet unnormahsed joint probabihty distribution of extensive quantities Q, P in the GCE is then 
given by the Fourier transform of Eq. (fT2|) after substitutions Eqs. (l21ll22l) : 



u) 



(2vr) 



exp 



V'^ ( (3, fl, u] 0, d 



(23) 



More details of the calculation, in particular on the connection between the partition functions 
Z'^'^iy, P, ll,u) and the conventional version ZMCE(y,Q,P) [42, |43|, can be found in Appendix Rl 
Depending on the system under consideration, we introduce the vector of extensive quantities Q and 
corresponding Wick rotated fugacities 0: 



Q = {B, S, Q,Na,Nb) 
Q = {Q,Na,Nb) 



B, <PS, (PQ-, (PNai (PNb, 



for a HRG , 
for a pion gas 



(24) 
(25) 



Here B is the baryon number, S is the strangeness, and Q is the electric charge of the system. 
Together with particle numbers Na and Nb this would be a 5-dimensional distribution in the case of 
a CE HRG. Additionally for four-momentum conservation, yielding a 9-dimensional Fourier transform 
Eq.(l23]) for a MCE HRG, we write: 



a 



{oiE, ap,, apy, apj 



(26) 



where E is the energy and Px, Py, and Pz are the components of the collective momentum of the 
system, while a are the corresponding fugacities. 



The integrand of Eq. fl2^ is sharply peaked at the origin = a = in the TL [4l|. The main 
contribution therefore comes from a very small region. To see this, a second derivative test can be 
done on the integrand of Eq. (!23l) taking into account the first two terms of Eq. (!29|) . The limits 
of integration can hence be extended to ±oo. The distinction between discrete (Kronecker 6) and 
continuous quantities (Dirac 6) is not relevant for the TL approximation, where particle number is 
a continuous variable to be integrated over. We thus proceed by Taylor expansion of Eq.( fT3l) . For 
this it is convenient to include everything into a common vector notation: 



Q 



Q, P 



and 



d 



0, d 



(27) 



The dimensionality of the vector Q is denoted asJ = 2 + 3 + 4 = 9fora MCE HRG. We now expand 
the cumulant generating function, \l/ ( /3, /I, m; 6' ) , in a Taylor series: 



E 

n=0 



w. 



9jl 9 32 ■ ■ ■ 9 in 1 



(28) 



where the elements of the cumulant tensor, k^^'-^^'' ■ are defined by: 



-I) 



de,,de,,...de,^ 



(29) 



Generally cumulants are tensors of dimension J and order n. The first cumulant is then a vector, 
while the second cumulant is a symmetric J x J matrix. A good approximation to Eq. (!23|) around 
the point Qeq = {Qeq, Peg), cau be found in terms of a Taylor expansion of Eq. (fT3l) in 6' = (0, a), if: 

dZ^{V,P,fl,u) 



dQ 



. 



(30) 



Implicitly, Eq.( l30i) does not define chemical potentials p, and four-temperature P = Pu, but 
corresponding Lagrange multipliers which maximise the amplitude of the Fourier spectrum of the 
generating function for a desired value of Qeq and Peq. Their values generally differ from the GCE 
set (/3, /i, , however they coincide in the TL. Lagrange multipliers can be used for finite volume 
corrections 411] • In the following we restrict ourselves to the large volume approximation. 



C. Joint Distributions 

In the large volume limit, i.e. V —>■ oo, one may use the asymptotic solution, and only consider 
the first two cumulants, Eq. (l29l) . Substituting Eq. (l28|) into Eq. (!23|) yields a standard J dimensional 
Gaussian integral with solution: 



ZQiV,T,fI,u) 



'GCE 



{V, (3, fl, u) 



deter 

where the elements of the new variable ^ are defined by: 

e = im' {^-% ■ 



exp 



(31) 



(32) 



The elements of the vector (AQ) measure the distance of the charge vector Q to the GCE mean 

{^Qt = {Q'-VkD , (33) 
and (cr^^) is the inverse square root of the second order cumulant K2- 

- ..-1/2 



(34) 



The GCE joint distribution of extensive quantities Q is a MND^: 

(35) 



'^'^^^ ZacEiV^P.n.n) - Cl^vyi^ deter 



Mean values in the TL are given by the first Taylor expansion terms, (A^a) = V . {Q) = Vk^, 
(E) = Vnf, etc. and converge to GCE values. To obtain a joint (two-dimensional) particle multi- 
plicity distribution one has to take a two-dimensional slice of the (J dimensional) GCE distribution, 
Eg. (135!) . around the peak of the extensive quantities which one is considering fixed. The co- variance 
tensor K2 will be spelled out exphcitly and discussed in Section lTlIBI Further details of the calculation 
can be found in Appendix [Bl 

III. FLUCTUATIONS AND CORRELATIONS WITHIN A MOMENTUM BIN 
A. Static System 

Let us start by discussing properties of a static thermal system. We want to measure joint 
distributions of multiplicities Na and Nb in limited bins of momentum space (of width ApT for 
transverse momentum bins and Ay for rapidity bins). Depending on the size and positions of the 
bins, one finds different fluctuations and correlations. Results will, in particular, be compared to 



the acceptance scaling approximation employed in 



31 



36|, 



33], which assumes random observation 

so Appendix [C]) . 



of particles with a certain probability q, regardless of particle momentum (see a 
Corresponding results for scaled variance in MB statistics can also be found in 29 1. 

For our examples we chose a gas with three degenerate massive particles (with positive, negative 
and zero charge) with mass m = 0.140GeV in three different statistics (MB, FD, BE). The momentum 
spectra are assumed to be ideal GCE spectra, due to the large volume approximation. In Fig.([T]) 
transverse momentum and rapidity spectra are shown for MB statistics. BE and FD statistics yield 
similar spectra, unless chemical potentials are large. 

We can then deflne bins by requiring each bin to hold the same fraction of the total multiplicity. 
Note that in this case the width and position of bins Ap^ and Ay will strongly depend on the 
underlying momentum spectra. Our examples, in particular the FD little academic in 

the sense that there is no fermion of this mass. In a HRG, often applied to heavy ion collisions, the 
lightest fermion is the nucleon for which quantum effects are probably negligible. 

In Fig. ([2]) we present the scaled variance u), calculated using Eq. lJTOl) . within different transverse 
momentum bins ApT (left) and rapidity bins Ay (right). The scaled variance in limited bins of 
momentum space is more sensitive to the choice of particle statistics than the spectra would suggest. 



Finite volume corrections to Eq. ([35| converge like in the TL 



m 




FIG. 1: (Color online) Differential particle spectra for a 'pion gas' at T = 160MeV. Transverse momentum 
spectrum (left) and rapidity spectrum (right). Both curves are normalised to unity. The bins are constructed 
such that each bin contains 1/9 of the total yield. 

BE and FD effects are particularly strong in momentum space bins in which occupation numbers are 
large. Hence, at the low momentum tail one finds suppression of fluctuations for FD and enhancement 
for BE, while at the high momentum tail, one finds ube — ^mb — ^fd, Fig.(l2]) (left). The rapidity 
dependence. Fig. ([2]) (right), has a different behaviour. The reason is, that in any Ay bin there 
is some contribution from a low pr tail of the differential momentum spectrum dN/dy/dpT where 
quantum statistics effects are pronounced. This leads to a clear separation of the curves and one 
finds ujbe > ^mb > ^fd- In contrast to this, the in integrated (all particles observed) scaled 



variance is rather insensitive to the choice of statistics |39l] (unless chemical potentials are large). 
Please note that there are in fact 3 different 'acceptance scaling' lines in Fig. ([2]), which extrapolate 
the An integrated scaled variance to limited acceptance. The differences are however very small and 
all 3 lines lie practically on top of each other. 

In Fig. ([3]) we present the correlation coefficient p, calculated using Eg. (ITT]) , between positively and 
negatively charged particles in transverse momentum bins Apx (left) and rapidity bins Ay (right). 
The 4n integrated correlation coefficient between positively and negatively charged particles would 
be /947r = 1 in the CE and MCE. In the GCE it would be 0. In the MB CE it would not show any 
momentum space dependence and would always be p > 0. In the MCE the situation is qualitatively 
different: in low momentum bins particles are positively correlated, while in high momentum bins 
they can even be anti-correlated. Horizontal lines again indicate acceptance scaling (Appendix [C|) . 
Quantum effects for the correlation coefficient remain small as there is no explicit local (quantum) 
correlations between particles of different charge. 

It should be stressed that the ApT dependence in Figs.( l2|3il is a direct consequence of energy 
conservation. The Ay dependence of u and p, however, is due to joint energy and longitudinal 
momentum (P^) conservation. Disregarding conservation leads to a substantially milder Ay 




P,[GeV] y 

FIG. 2: (Color online) Transverse momentum (left) and rapidity dependence (right) of the MCE scaled 
variance of negatively charged particles at T = 160MeV, for a MB (blue), FD (green), BE (red) 'pion gas' 
at zero charge density. Momentum bins are constructed such that each bin contains the same fraction q 
of the average tt~ yield. The horizontal bars indicate the width of the ApT or Ay bins, while the marker 
indicates the position of the center of gravity of the corresponding bin. Dashed lines indicate acceptance 
scaling results, Eq. (|C6p . 
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FIG. 3: (Color online) Same as Fig. ([2]), but for the MCE correlation coefficient between positively and 
negatively charged particles. Dashed lines indicate acceptance scaling results, Eq. (|C9|) . 



dependence, see Fig.(jl]). 

This behaviour can be intuitively explained: in a low momentum bin it is comparatively easy to 
balance charge, as each individual particle carries little energy and momentum. In contrast to this, 
in a high momentum bin with, say an excess of positively charged particles, it is unfavourable to 
balance charge, as one would also have to have more than on average negatively charged particles, 
and each particle carries large energy and momentum. This leads to suppressed fluctuations and 
correlations in high momentum bins when compared to low momentum bins. 

In a small mid-rapidity window, with \y\ < 0.3, the effects of globally applied motional conserva- 
tion laws cease to be important (see Fig. ([5])). Local correlations due to BE and FD statistics begin 




FIG. 4: (Color online) Rapidity dependence of the scaled variance of positively charged particles (left) and 
of the correlation coefficient between positively and negatively charged particles (right). Calculations are 
done for the same system as in Figs. (|2l3p . however disregarding momentum conservation. 
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FIG. 5: (Color online) Transverse momentum dependence of MCE scaled variance of negatively charged 
particles {left) and the MCE correlation coefficient between positively and negatively charged particles 
{right). Only particles in a mid-rapidity window \y\ < 0.3 are measured. Dashed lines denote the GCE 
result, Eq.dM])- 



to dominate, and MCE deviations from the GCE results, Eq. (!36l) . are relevant only for the highest 
momentum bins. In BE or FD statistics we find for vanishing bin size {6y, 6pt)'- 



GCE 
^5 



Ns,Ns 



1 



-|- g—/3mT cosh 



(36) 



BE and FD effects are strongest around mid-rapidity y = 0. MCE calculations in Fig. ([5]) are close 
to the GCE estimate Eq. fl36|) . In MB statistics we find only a weak Ap-p dependence in a small mid- 
rapidity window. Please note that the acceptance scaling procedure predicts a Poisson distribution 
with uj ^ 1 and p ^ for all three statistics in the limit of very small acceptance. 



B. Collectively Moving System 



As established a long time ago, in order to properly define the thermodynamics of a system with 



collective motion, the partition function needs to be Lorentz invariant |48l . |49j. The expectation 
values of observables need hence to transform according to the Lorentz transformation properties 
of these observables. In particular, the temperature T is promoted to a four-vector jS'^ (combining 
local temperature with collective velocity). The entropy, as well as particle multiplicities, remain 
Lorentz-scalars. 

These requirements are in general not satisfied unless momentum conservation is put on an equal 
status with energy conservation. If the system is described by a MCE, then momentum should be 
conserved as well as energy 48|, |49l]. If the system is exchanging energy with a bath, the bath needs 
to exchange momentum as well. 

For ensemble averages, neglecting these rules and treating momentum differently from energy 
is safe as long as the system is close to the thermodynamic limit, since there ensembles become 
equivalent. The same is not true for fluctuation and correlation observables, which remain ensemble- 



specific 



For a system at rest, these requirements are not apparent since the net momentum is zero. Statis- 
tical mechanics observables in a collectively moving system, however, lose their Lorentz invariance, 
if this is not maintained in the definition of the partition function. 

To illustrate this point, we consider the properties of a system moving along the 2;-axis with a 
collective velocity given by Eq. ffTSl) . The total energy of the fireball is then E = Mcosh(?/o), while its 
total momentum is given hj = M sinh(?/o)- The mass of the fireball in its rest frame is M = P^u^. 
The system 4-temperature is (3 = f3u. Local temperature and chemical potentials remain unchanged. 
We will use this section for a discussion of the second rank tensor (or co- variance matrix) ^2, Eq.f l29p . 

The second order cumulant ^2, Eq.f l29l) . is given by the second derivatives of the cumulant generat- 
ing function with respect to the fugacities. Essentially this is the Hessian matrix [s^ of the function 
Eq.f ll3l) . encoding the structure of its minima. The diagonal elements the variances of the 

GCE distributions of extensive quantities X. For example, k^'^'^'^ measures the GCE variance of 
the distribution of particle multiplicity of species A, while k^'*^ denotes the GCE electric charge 
fluctuations, etc. The off-diagonal ^2 ' elements give GCE co-variances of two extensive quantities 
X and Y. 



For a boost along the z-axis the general co-variance matrix for a pion gas reads: 
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(37) 



Off-diagonal elements correlating a globally conserved charge with one of the momenta, i.e. k^'^"", 

p p 

as well as elements denoting correlations between different momenta, equal to zero 

due to antisymmetric momentum integrals. The values of elements correlating particle multiplicity 
and momenta, i.e. k^^'^"" ^ depend strongly on the acceptance cuts applied. For fully phase space 
integrated (47r) multiplicity fluctuations and correlations these elements are equal to 0, again due to 
antisymmetric momentum integrals. 

It is instructive to review the transformation properties of k under the Lorentz group: Kg ' 
contains the correlations between 4-momenta and, in general, (scalar) conserved quantities and 
particle multiplicities . Hence, the elements k^'''^ , i.e. (AP^AP'^), will transform as a tensor 
of rank 2 under Lorentz transformations; K2^'^\ i.e. (AP^AQ^), will transform as a vector (the 



rapidity distribution will simply shift); and the remaining Kg will be scalars 



For a static system one finds for the co-variances k 



E,Px 



E,P, 



E,P, 



0. Under these 



two conditions, a static system and full particle acceptance, the eigenvalues of the matrix Eq. (lH7I) 
factorise, and momentum conservation can be shown to drop out of the calculation 29 1. 

For a boost along the z-axis (and arbitrary particle acceptance) it is the appearance of non- 



vanishing elements 

^2 — ^2 

' 7^ which make the determinant of the matrix Kg, Eq. (!37|) . invariant 



E,P, 



against such a boost. Please note that still k 



Px,E 



Kj'-) 



0. 



In Fig. ([6]) we show multiplicity fluctuations (left) and correlations {right) for a system with boost 
Uq = 2. The rapidity spectrum of Fig.([l]) (right) is simply shifted to the right by two units. The 
construction of the acceptance bins is done as before. Multiplicity fluctuations and correlations 
within a bin transform as a vector (i.e., its z component shifts in rapidity) as inferred from their 
Lorentz-transformation properties, provided both energy and momentum along the boost direction 
are conserved. 

I — I 

This last point deserves attention because usually (starting from [17]) micro canonical calculations 
only conserve energy and not momentum. Imposing exact conservation for energy, and only average 
conservation of momentum will make the system non-Lorentz invariant, since in a different frame 
from the co-moving one energy and momentum will mix, resulting in micro state-by-micro state 




FIG. 6: (Color online) Same as Fig. ([2]) (right) and Fig. ([3]) (right), but for a system moving with yo = 2. 



fluctuations in both momentum and energy^ 




FIG. 7: (Color online) Same as Fig. ([6]), but without Pz exact conservation. 

This situation is explicitly shown in Fig. ([7]). Here, we have calculated multiplicity fluctuations 
and correlations in the same system as in Fig. ([6]), but with exact conservation of only energy (and 
charge). In the co- moving frame of the system, the fluctuations and correlations are identical to 
Fig.(jl]). When the system is boosted, however, the distribution changes (not only by a shift in 

^ This result is somewhat confusing, because energy-momentuni is a vector of separately conserved currents. It 
is therefore natural to assume that these currents can be treated within different ensembles; they are, after all, 
conserved separately. It must be kept in mind, however, that it is not energy-momentum, but particles that are 
exchanged between the system and any canonical or grand canonical bath. The amounts of energy and momentum 
carried by each particle are strictly correlated by the dispersion relation [49]. In the situation examined here 
(unlike in a Cooper- Frye formalism [si] ]. where the system is "frozen" at the Freeze-out hypersurface, a space- 
time 4-vector correlated with 4-momentum) all time dependence within the system under consideration is absent 
due to the equilibrium assumption. Furthermore, the system is entirely thermal: the correlation between particle 
numbers when the system is sampled "at different times" is a (5— function, that stays a (5— function under all Lorentz- 
transformations. Hence, unlike what happens in a Cooper-Frye freeze-out, energy-momentum and space-time do 
not mix in the partition function. Together with the constraint from the particle dispersion relations, this means 
that different components of the 4-momentum need to be treated by the same ensemble, as explicitly demonstrated 
in this section. 



rapidity, as required by Lorentz-invariance), and loses its symmetry around the system's average 
boost. 

This last effect can be understood from the fact that momentum does not have to be conserved 
event by event, but energy does. It is easier, therefore, to create a particle with less rapidity than 
average (having less momentum than the boost, and parametrically less energy) than with more 
rapidity than average (having more momentum than the boost, and parametrically more energy) 
and still conserve energy overall. This leads to suppressed multiplicity fluctuations and a negative 
correlation coefficient for rapidity bins in the forward direction in comparison to rapidity bins in 
the backward direction. In Fig. ([6]), where the system also needs to conserve momentum exactly, 
this enhancement is balanced by the fact that it will be more difficult to conserve momentum when 
particles having less momentum than the boost are created. 

A situation such as that in Fig.Q is impossible to be realized physically. It could, however, 
be realized within "system in a box" -type calculations with non-equilibrium models: E.g., a trans- 
port model inside an infinitely heavy box (that absorbs momentum but not energy event-by-event) 
would end up exhibiting micro canonical correlations similar to those in Fig.([7j). A similar box with 
'periodic' walls, however, would conserve energy as well as momentum inside the box, and should 
therefore behave as in Fig.(l6l). Thus, correlations within boosted sources provide a sensitive test of 
the Lorentz-invariance of such transport models. 



IV. CORRELATIONS BETWEEN BINS DISCONNECTED IN MOMENTUM SPACE 



"Long range correlations" between bins well disconnected in momentum have been suggested to 
arise from dynamical processes. Examples include color glass condensate [?, 0], droplet formation 



driven hadronization [s^, and phase transitions within a percolation-type mechanism 53. 54 



elliptic flow measurements, widely believed to signify the production of a liquid at RHIC 



55 



'he 



57 



56, 



are also, ultimately, correlations between particles disconnected in phase space (here, the 
azimuthal angle). 

As we will show, however, conservation laws will also introduce such correlations between any 
two (connected or not) distinct regions of momentum space. No dynamical effects are taken into 
consideration (only an isotropic thermal system). 

Let us first consider correlations between the multiplicities of particles A and B, within two bins 
centered around i/a and i/b, with (constant) widths Ai/a = Ai/b = 0.2. In Fig.® (left) we show the 
correlation coefficient, calculated using Eq.l ITT]) . between positive and negative particles as a function 
of HA and i/b- In Fig.® [right) we show the correlation coefficient between like-charge, unlike-charge, 
and all charged particles as a function of i/gap = Va — Vb- 

Energy conservation always leads to anti-correlation between different momentum space bins. 




A 'gap 



FIG. 8: (Color online) Left: The correlation coefficient p between the multiplicity of positively charged 
particles in a bin located at yA with negatively charged particles in a bin located at i/b, both with a 0.2 
width in rapidity. Right: The correlation coefficient between multiplicities in two bins separated by ygap of 
like, unlike, and all charged particles. Both plots show MCE MB results. 

Charge conservation leads to a positive correlation of unlike charged particles and anti-correlation of 
like-sign particles. Longitudinal momentum conservation, however, is responsible for the structure 
visible in Fig. (left). Having a small (large) number of particles in a bin with positive average 
longitudinal momentum, leads to a larger (smaller) number of particles in a bin with different but 
also positive P^, (blue dips). This makes also a state with smaller (larger) particle number with 
opposite average longitudinal momentum —Pz more likely (red hills). At large values of t/A the 
correlation coefficient p ~ for any ys, because the yield (Na) in Ay a is asymptotically vanishing. 

In Fig.(l8]) (right) we show the correlation coefficient along the diagonal from top left to bottom 
right as a function of separation. Unlike-sign particles are positively correlated. Like-sign and 
all charged particles are negatively correlated at small separation ygap. For large separation the 
correlation becomes asymptotically zero, because the yield is zero. However, please note that in 
particular p{n'^,Tr'^) > p{Tr~^,n~) at large ygap. Pz conservation is dominant. 

Disregarding Pz conservation would destroy the particular structure in Fig. ([8]) (left) and lead to 
a single peak at the origin. The correlation would then be insensitive to the momentum direction, 
and only be sensitive to the energy content of a bin Ay. The observables in Fig. ([8]) transform under 
boosts {yA,B yA,B — Ho), provided momentum along the boost axis is exactly conserved. 

Angular correlations also arise due to conservation of transverse momenta Px and Py. In Fig. ([9]) we 
show the correlation coefficient between particles in different A0 bins. The flat^ angular spectrum 
dN/d(j) has been divided into 10 equal size bins and the correlation coefficient is presented as a 
function of separation of the centers of the corresponding bins. 



^ Since we consider globally equilibrated systems, elliptic flow is disregarded here. 
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FIG. 9: (Color online) The correlation coefficients of particles in distinct A(j) bins as a function of separation 
(pgap in azimuth, {left) integrated over all phase space, {right) only particles with \y\ < 0.3 are observed. 
Both plots show MCE MB results. No elliptic flow is considered. 

To explain Fig.dH]) we note firstly that when disregarding exact conservation of Px and Py the 
correlation coefficients are insensitive to the distance (pgap of any two bins. Only the correlations 
due to energy and charge conservation affect the result. Charge conservation leads to correlation of 
unlike-sign particles and to anti-correlation of like-sign particles. Energy conservation always anti- 
correlates multiplicities in two bins. For p('7r^,7r^) the effect of charge conservation cancels for a 
neutral system, however, effects of energy-momentum conservation are stronger, as a larger number 
of particles (hence a larger part of the total system) is observed. 

Conservation of transverse momenta P^ and Py is now responsible for the (pgap dependence of p. 
The line of arguments is similar to the ones before: Observing a larger (smaller) number of particles 
in some bin at 0o implies that, in order to balance momenta P^ = Py = 0, one should also observe 
a larger (smaller) number of particles in the opposite direction vr — 0o- A larger (smaller) number of 
particles in a bin with (pgap = tt/2 would do little to help to balance momentum, but conflict with 
energy conservation. 

V. DISCUSSION AND SUMMARY 

We have presented multiplicity fluctuations and correlations in limited momentum bins for ideal 
relativistic gases in the MCE in the thermodynamic limit. For our examples we chose a gas with 
three degenerate massive particles (positive, negative, neutral) in three different statistics (Maxwell- 
Boltzmann, Fermi-Dirac, Bose-Einstein). 

For the width of multiplicity distributions in limited bins of momentum space a simple and 
intuitive picture emerges. In the Maxwell-Boltzmann approximation one finds a wider distribution for 
momentum bins with low average momentum when compared to bins with higher average momentum, 
but same average particle number. This qualitative behaviour is a direct consequence of energy and 



momentum conservation. The results in Fermi-Dirac and Bose-Einstein statistics, furthermore, show 
pronounced effects at the low momentum tail of the spectrum. 

The correlation coefficient additionally shows a similar qualitative behaviour. In bins with low 
average momentum the correlation coefficient between positively and negatively charged particles 
is indeed positive, as one would expect from charge conservation. However, in bins with large 
average momentum the effects of joint energy and momentum conservation can lead to anti-correlated 
distributions of unlike-charged particles. 

For boosted systems we found that the role of exactly imposed motional conservation laws is 
particularly important. Fluctuations and correlations transform under boosts, provided momentum 
conservation along the boost direction is taken into account. This ensures, in particular, that they 
become boost invariant if the underlying system is boost-invariant. 

Lastly, we found that even in the thermodynamic limit long range correlations between discon- 
nected regions in momentum space prevail. Multiplicities in different rapidity bins, as well as different 
bins in azimuth, have a non-zero correlation coefficient. 

It is premature to use the model presented here for a quantitative comparison with experimental 
data. Firstly, the inclusion of resonances will provide important corrections. These can be imple- 
mented within our model using Monte-Carlo techniques, and will be the subject of a subsequent 



work 



59|. 



An additional effect missing here that could change results qualitatively is longitudinal flow. It 
can be seen, "by symmetry", that correlations and fluctuations in a perfectly boost-invariant fluid 
would be, just like other physical observables, independent of rapidity. Reconciling this with the 
calculations in sections IIII Al and IIIIBI would require calculating correlations of many independent 
sources, each centered around a particular rapidity. 

It is however not currently clear whether experimental data, even at highest RHIC and LHC 
energies at mid-rapidity, approximates this limit. At lower SPS energies, measurements [id] of the 
rapidity and transverse momentum dependence of particle multiplicity fluctuations show qualitatively 
similar results to those of our calculations (it should be noted that, as shown in |6Ci], this behaviour for 
fluctuations also arises in molecular dynamics models, where conservation laws are included but equi- 
librium is not assumed). It has long been noticed 6ll, 1621] that many observables binned in rapidity 
obey "universal fragmentation" , suggesting that a "Landau hydrodynamics" , with negligible initial 
longitudinal flow even away from mid-rapidity, might be more appropriate than boost-invariance 
to describe the initial state, even at ultra- relativistic energies |63|. Experimentally, these measure- 
ments might be non-trivial since the size of the bins must be lar ge e nough for conservation laws 
to have an effect, but, since RHIC has large- acceptance data 55|, l56| and the LHC is planning a 
larger-acceptance detector js^, they are in principle possible. 



If boost-invariance is not really there, then correlations and fluctuations binned by rapidity should 
be qualitatively similar to those calculated in section IIII Al even at RHIC and LHC energies. We 
therefore suggest the experimental measurement of such fluctuations (Figs. (!2|3|) (right)) at these 
energies as an experimental probe of the degree of boost-invariance of the system. 

Similarly, transverse flow should not change the correlation and fluctuations within limited trans- 
verse momentum bins (the left panels of Figs. fl2][51) ) beyond a trivial shift provided there are no 



significant fiuctuations in co^ 



lective fiow observables. These fiuctuations are widely expected to 



arise ■„ an .mperfect fluid M, but have remarkably uot beeu observed, for example, m ellipfc flow 

measurements [66|, |67|]. We suggest, therefore, that a measurement of pt binned fiuctuations and cor- 
relations could provide a qualitative way to assess the magnitude of event-by-event fiow observables 
wrt to the thermal observables presented here. 

Finally, effects of the sort studied in section IIVI will surely appear in any measurement looking 
for correlations across momentum space. Long range rapidity correlations have, in particular, been 
advocated as a signature of new physics l7, 81. Our work shows that correlations due to conservation 

nnn 

laws actually cover as wide a range in rapidity as those measured in [12, [13!, [IJ] . The magnitude of 
the correlations, however, is significantly (as much as an order of magnitude) lower than either the 
experimental result or any reasonable "new physics" , suggesting that the origin of the experimentally 
observed correlations lies elsewhere. In particular, correlations induced by initial state geometry are 
considerably larger than those induced by conservation laws, as a comparison between Fig.® and 
the results in [151] will show^. Nevertheless, energy-momentum conservation does trigger correlations 



across a wide rapidity interval, and, as shown in 



60|, qualitatively the magnitude of these correla- 



tions is independent of the degree of equilibration of the system, so their presence in experimental 
data is very plausible. Perhaps complementing the rapidity correlations with azimuthal correlation 
measurements, such as those in Fig.Q, might clarify their role, although the latter are particularly 
susceptible to "non-trivial" physics contributions, such as jet pairs and elliptic fiow. 

In conclusion, we have presented microcanonical ensemble calculations of correlation and fiuctu- 
ation observables within and across bins within a range of rapidity and transverse momenta. The 
calculations presented here provide qualitative effects affecting multiplicity fluctuations and correla- 
tions. These effects arise solely from statistical mechanics and conservation laws. It will be extremely 
important to see whether these qualitative effects are visible in further experimental measurements 
of the momentum dependence of multiplicity fluctuations and correlations. If so, these effects might 
well be of similar magnitude to the signals for new physics. Disentangling them from dynamical 



^ Note, however, that these results were obtained in the infinite volume limit. Finite volume effects are likely to increase 
the strength of correlations arising from conservation laws, though for realistic nuclear volumes such corrections 
should not alter the results by an order of magnitude 



correlations will then be an important, and likely non-trivial task. 
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APPENDIX A: MCE PARTITION FUNCTION 



This section serves to provide a connection between Eqs.(l3]l5l) and Eq. (l23i) : namely to prove the 
following relation: 

Z«'^^(V^, (3, fl, u) = e-^""'^^ e'i'^^P Zmce{V, P, Q) , (Al) 

where Zmce{V, P,Q) is the standard MCE partition function for a system of volume V, collective 
four- momentum P and a set of conserved Abelian charges Q, as worked out in 42, |43l]. The MCE 
partition function Zmce{V, P, Q) counts the number of micro states consistent with this set of fixed 
extensive quantities. Likewise one could interpret the number Z'^'^CV, (3, fl,u) as the number of 
micro states with the same set of extensive quantities for a GCE with local inverse temperature (3, 
four-velocity u, and chemical potentials /I. 

The starting point for this calculation is our Eq. fl23l) : 



(2^ 



exp 



V'^ ( /3, /2, u] 0, a 



(A2) 



Let us take a closer look at the exponential of Eq. flA2l) . For this we spell out Eq.f llSp and use the 
substitutions Eqs.(l2ll) and ([221): 



exp 



±1 



(A3) 



Expanding the logarithm yields: 

Vgi 



exp 



E 



ni 



(A4) 



Replacing now the momentum integration in Eg. (]A4|) by the usual summation over individual mo- 
mentum levels f cPp — > J2kn Si'^^s: 



exp 



^ ^ ^ rii *^ 



(A5) 



Finally, expanding the exponential yields: 



(2^ 
1 f 91 iTl) 
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nnnE 



Only sets of numbers {cfc„j } which meet the requirements: 



(A6) 



(A7) 



E E E = E E E = 

I ni=l knj I ni=l 

have a non-vanishing contribution to the integrals. Therefore we can pull these factors in front of 
the integral: 

■K OO 



u 



oo oo 



nnnE 



1 fgdTi) 
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I ni=l Cfe„^=0 "-"i 
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rii 



Reverting the above expansions one returns to the definition of ZMCEiV, P,Q) from Ref.[42|, |43 | 
times the Boltzmann factors: 



d^a 



exp 



{2ny J {2nr 

— TV —OO 



which proves Eq. llAip . Therefore we write for the GCE distribution of extensive quantities: 

^ ZGCE{V,P,fI,u) Zgce{V,P,11,u) ' 

which provides the promised connection between Eqs.([3]l5]) and Eq. (l23l) . 
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APPENDIX B: ASYMPTOTIC JOINT DISTRIBUTION 



The MCE joint multiplicity distribution Pmce{NA, Nb) is conveniently expressed by the ratio of 
two GCE joint distributions: 

(Bl) 
(B2) 



PmceiNA,NB) = P,UNa,Nb\B,S,Q,E,...) 

_ P,UNa,Nb,B,S,Q,E,... 



PgceiB, S,Q, E, . . .) 



In the TL the distributions Pgce{NA, Nb, B, S,Q, E, . . .) and Pgce{B, S,Q, E, . . .) can be approxi- 
mated by MND's, Eg. (13 II) . The charge vector Eg. (1331) for a MCE HRG with three conserved charges 
would read: 

(AQ) = (ANa, ANb, AB, as, AQ, AE, ...) • (B3) 
Evaluating the MND, Eg. (!3T]) . around its peak for {B, S, Q, E, . . .) yields: 



(AQ) = {ANa, ANb, 0, 0, 0, 0, 
The vector Eg. (|32ll then becomes: 

/ Ai,i AA^^ + Ai,2 ANb\ 
A2,i AA^^ + A2,2 ANb 
A3,i AA^^ + A3,2 ANb 
A4,i AA^^ + A4,2 ANb 
A5,i ANa + A5,2 ANb 



(B4) 



\ 

where Xij are the elements of the matrix Eg. fl34|) . Therefore: 



(B5) 



J J J 

= [{ANAfJ2 4i + 2 (AA^^) (AAT^) ^ A,,iA,,2 + (AAr^)'^^ 
j=i j=i j=i 



(B6) 



with J = 2 + 3 + 4 = 9 for a MCE HRG with momentum conservation. Using Eg.( lB6p . the micro 
canonical joint multiplicity distribution of particle species A and B can thus be written as: 



PmceiNA,NB) 



1 detdAT 



exp 



2 



(B7) 



{2nV) det a 

where is the 7-dimensional inverse sigma tensor of the distribution Pgce{B, S,Q, E, . . .). Com- 
paring this to a bivariate normal distribution Eg. , one finds: 
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After a short calculation one finds for the co-variances: 
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(Bll) 
(B12) 
(B13) 



and additionally the correlation coefficient, Eq. (fTT|) : 

p = ^ = ^ , (B14) 

where the terms A,B,C are given by Eqs. (]B8l - IBlOp . For the normalization in Eq. (lB7|) (from a 
comparison with Eq.([6])) one finds: 

det (Ttv 1 



deta aAasVil-p^) 
APPENDIX C: ACCEPTANCE SCALING 



VAB . (B15) 



To illustrate the 'acceptance scaling' procedure employed in [3l|, l36|, l37|] we assume uncorrected 
acceptance of particles of species A and B. Particles are measured or observed with probability q 
regardless of their momentum. The distribution of measured particles ua, when a total number Na 
is produced, is then given by a binomial distribution: 

PacciriAlNA) = g"-^ (1 - ^)'^""''" (r^) • (CI) 

The same acceptance distribution is used for particles of species B. Independent of the original 
multiplicity distribution P{Na, Nb), we define the moments of the measured particle multiplicity: 

K-n^) = Yl Yl n'Xn'^PaccinAlNA) PacciriBlNB) P{Na,Nb) . (C2) 

nA,nB Na,Nb 

For the ffist moment (ua) one finds: 

{ua) = q {Na) • (C3) 
The second moment {n\) and the correlator {ua ■ ub) are given by: 

{n\) = q'{Nl) + q{l-q){NA) , (C4) 
{ua-ub) = q^NA-NB) . (C5) 

For the scaled variance uj^ of observed particles one now finds [311] : 



~ W) - ~ 9 + 9 ^47r ) I'-'OJ 



where a;^ is the scaled variance of the distribution if all particles of species A are observed. Lastly, 
the correlation coefficient pq is: 

Pq = I . 

^J{{/\uAf) {(Aub)') 

with {AuaAub) = {ua ■ ub) — {ua){ub), and {{Aua)"^) = {u\) — {ua)^- Substituting the above 

relations, one finds after a short calculation: 

/ r 1 ^1/2 

Pg = Pa. q \ <^l + g(i - q)^t + g(i - qX + (1 - qf ■ (C8) 



In case ojf^ = uf^ = u^.,,, Eq. OCSP simplifies to: 



q + q uJ4n 

Both lines are independent of the mean values (Na) and (Nb). 



Pq = PiTT Z ■ . C9 

1 - q + q ojA-jr 
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